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We examine transient axial creeping flow in the annular gap between a rigid cylinder 
and a concentric elastic tube. The gap is initially filled with a thin fluid layer. The study 
focuses on viscous-elastic time-scales for which the rate of solid deformation is of the same 
order-of-magnitude as the velocity of the fluid. We employ an elastic shell model and the 
lubrication approximation to obtain a forced nonlinear diffusion equation governing the 
viscous-elastic interaction. In the case of an advancing liquid front into a configuration 
with a negligible film layer (compared with the radial deformation of the elastic tube), 
the governing equation degenerates into a forced porous medium equation, for which 
several closed-form solutions are presented. In the case where the initial film layer is non- 
negligible, self-similarity is used to devise propagation laws for a pressure driven liquid 
front. When advancing external forces are applied on the tube, the formation of dipole 
structures is shown to dominate the initial stages of the induced flow and deformation 
regimes. These are variants of the dipole solution of the porous medium equation. Finally, 
since the rate of pressure propagation decreases with the height of the liquid film, we show 
that isolated moving deformation patterns can be created and superimposed to generate 
a moving wave-like deformation field. The presented interaction between viscosity and 
elasticity may be applied to fields such as soft-robotics and micro-scale or larger swimmers 
by allowing for the time-dependent control of an axisymmetric compliant boundary. 


1. Introduction 

We examine the effect of elasticity on transient axial creeping flow in the annular gap 
between a rigid cylinder and a concentric elastic tube. Specifically, we aim to analyze 
nonlinear viscous-elastic dynamics where the elastic deformation significantly changes 
the boundaries of the flow. The gap between the cylinder and tube is assumed small 
compared to the radius and length of the cylinder. The flow-field is modeled by applying 
the lubrication approximation, while the external elastic tube is modeled by thin shell 
theory. We focus on viscous-elastic time-scales for which the rate of solid deformation is 
of the same order-of-magnitude as the velocit y of the flow . _ 

A similar configuration was examined by IPaidous.^ ( 19981) who focused on high 
Reynolds numbers. The interaction of elas tic tubes and thin viscous flo ws, at the limit of 
low Re ynolds numbers, has been studied bv lHaloern fc Grotberd ( 1992 1 and White fc Heill 
(I2n05h who, among others, analyzed the dynamics of a liquid film coating the inner sur¬ 
face of an elastic tube in the context of flows in small airways in the lungs. Other studies 
involving a thin creeping flow between a rigid surface and an elastic surface focused on 
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planar configurations. These include Chauhan fc Badk3 ( 2002 1. who modeled the dy¬ 
namics of a contact lens during bl inking as a thin viscous film co ntained between an 
elastic shell and a flat rigid surface. Pihler-Puzovic et~^ ( 2012 ) and lAl-Housseinv et al 


(|2013r) studied the effect of elasticity on t he on set of Taylor-Saffman fingering instability 
in Hele-Shaw cells. IPihler-Puzovic et al. ( 2014 ) rel ated the patterns of visco us fingering 
to patt erns of wrinkling in an elastic Hele-Shaw cell. lTrinh et all ( 2014cJ) and Trinh et all 
(12014611 studied rigid and elastic plate s, either pinned or fre e floating, moving over a vis¬ 
cous film laying on a flat rigid surface. Carlson et (1201,^ studied the deformation and 
flow-field created by a propagating adhesion front attaching an elastic sheet to a rigid 
surface, where the gap between the elastic sheet and rigid surface is filled with a viscous 
fluid. 

For the case in which there is initially a vanishingly small gap between the rigid cylinder 
and the elastic t ube, the examined problem can be viewed as a cylindrical version of th e 
peeling problem ( McEwan fc Tavloi 19661 : iHosoi fc Mahadevan 2004 ; Lister et '^ 120131) . 
In this regard the peeling for mation of the current study bears mathematical pro ximity to 
a viscous gravity current fe.g.. BuckmasteiJl977 : HuPDertlll982 : Momoniatll2006l) . Gravity 
currents have also been studied in geophysical contexts when c oupled with elastic surfaces 
( Howell et al. 20131 : iHewitt et ai]l2015t [Balmforth et al. 2015 ). 

To the best of our knowledge, the problem examined in this work has not previously 
been studie d. The suggested configuration may have bearing on models of compliant 
boundaries ( Gad-el Hakll2003) . axisymmetric swimmers and soft robotic applications. 
The structure of this work is as follows: In §2 the geometry, relevant parameters and phys¬ 
ical assumptions are defined. In §3 elastic shell theory and the lubrication approximation 
are employed to obtain a forced nonlinear diffusion equation governing the viscous-elastic 
interaction. §4 examines an advancing liquid front in the limit of an initial vanishingly 
small gap (viscous peeling) where the governing equation simplifies to a porous medium 
equation. In §5 the self-similar scheme of the porous medium equation is modified to 
account for a pre-wetted gap between the cylinder and tube. In §6 the dipole solution of 
the porous medium equation is related to the pressure and deformation field created in 
the initial stages of the response to an advancing external pressure acting on the elas¬ 
tic tube. §7 examines isolated moving deformations for the case of relative axial speed 
between the rigid cylinder and the elastic tube. §8 presents concluding remarks. 


2. Problem formulation 

We study Newtonian, incompressible axial creeping flow in the annular gap between 
a rigid cylinder and a concentric linearly elastic tube. The gap between the cylinder 
and tube is assumed small compared to the radius of the cylinder. The flow is driven by 
time-varying pressures at the inlet or outlet, relative axial speed between the centre-body 
(cylinder) and elastic tube, as well as an external pressure field acting on the elastic tube. 
We examine transient dynamics in the viscous-elastic time-scale, where the rate of elastic 
deformation is of the same order-of-magnitude as the fluid velocity. We assume sufficient 
height of the liquid filled gap between the tube and centre-body, so that Van-der-Waals 
forces between the solids can be neglected. While §4-7 deal with configurations with 
constant radii of the centre-body and tube when at rest, in §2-3 we allow for small axial 
variations of the configuration (see Fig. [T|). 

Hereafter, normalized variables are denoted by uppercase letters and characteristic 
parameters are denoted by lowercase letters with asterisks (e.g., if a is a dimensional 
variable, a* is the characteristic value of a and A = a/a* is the corresponding normalized 
variable). 
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Figure 1: Schematic illustration of the elastic tube and annular fluid layer enclosing 
a weakly varying cross-section rigid centre-body. We define the cylindrical coordinates 
(r, z). Pe is the external pressure field, and dr, dz are radial and axial solid deformations, 
respectively. rh{z) is the radius of the centre-body, r^o is the minimal value of rb{z) and 
Sb(z) = rb(z) — Tbo- ri(z) + dr(z,t) and ro(z) + dr(z,t) are the inner and outer surfaces 
of the elastic tube, respectively. ho(z) + dr(z,t) is the annular gap between the inner 
centre-body and the elastic tube. po(t) is an arbitrary time varying inlet pressure. 


The relevant variables and parameters are time t, axial coordinate z, radial coordinate 
r, axial liquid speed Uz, radial liquid speed Ur, liquid pressure p, liquid viscosity p., liquid 
density p, solid axial deformation dz, solid radial deformation dr, solid strain and stress 
(Tij (acting on the plane normal to coordinate i and in the direction of coordinate j), 
solid Young’s modulus E, solid Poisson’s ratio iz, external shell (tube) inner contour ri{z) 
and outer contour ro{z) (see Fig.[T]), external shell midsection, rmiz) = {ri{z) + ro{z))/2, 
centre-body contour radius rb{z) and the length of the configuration in the z direction, 
1 . 

The gap (at rest) between the inner body and tube is defined as ho{z) = ri{z) — rb{z). 
We define the auxiliary radial coordinate, s = r — rb{z). The base radius of the inner-body 
is defined by rbo = min(rb( 2 ;)) and the corresponding function describing the centre- 
body radius is Sb{z) = rb{z) — rbo- We define the slenderness ratio of the configuration as 
e = rbojl. The ratio of characteristic fluid layer thickness at rest to characteristic radial 
deformation is defined by \h = h^/d* and the ratio of mean centre-body depth in excess 
of the base radius to characteristic radial deformation is defined by Xb = sl/d*. 

We define normalized coordinates, {R, Z) = {r/rbo,z/l), normalized radial and axial 
deflections, {Dr, Dz) = {dr jd*r,dzld*^, normalized liquid pressure, P = p/p* and external 
domain pressure, Pe = Pe/p*, where p* is the characteristic pressure, and normalized 
liquid velocity {Ur,Uz) = {ur/u/,Uz/u*^), with corresponding characteristic values m* 
and u*. 

Small parameters in the analysis include the ratio between the change in centre-body 
cross-section and the length of the geometry in the 2 ; direction, 


^ ^ ^ « 1 
dz I ' 

and the ratio between the characteristic radial deformation d) 
geometry in the z direction, i.e., the fluidic aspect ratio. 


( 2 . 1 ) 

and the length of the 
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El <C 1 . 


( 2 . 2 ) 
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We utilize a thin elastic shell approximation and thus require the ratio of wall thickness 
ro — ri to inner tube radius to be small, 

= £2 « 1, (2.3) 


as well as small elastic deformations, expressed by the ratios, 

d* d* 

(2.4) 

Ti I 

Our analysis also assumes negligible inertia, represented by the condition, 

= (2.5) 

A more detailed discussion of the condition of negligible inertia is presented in §3.4. 


3. Obtaining the leading order governing equation of the 
viscous-elastic dynamics 

In order to obtain the governing viscous-elastic interaction equation, we apply the 
lubrication approximation for the flow field and the Kirchhoff-Love thin shell approxima¬ 
tion for the solid deformation field, under the requirement of similar time-scales of both 
elastic and fluidic dynamics. 


Wh ile the elastic shell configuration is similar to our previous work (jElbaz fc Gat 
2014h . the characteristic shear stress applied by t he liquid on the solid is 0 ( 63 ) smaller 
compared with the characteristic shear obtained in lElbaz fc Gat ( 2014h . (this is readily 
seen from order-of-magnitude analysis of an integral force balance p*d*2'Kri ^ a*j.l2'Kri). 
We thus redevelop the thin shell approximation herein. 

3.1. The elastic problem 

The deformation field of axisymmetric linearly elastic material with negligible inertia is 
governed by the momentum equations, 


d d 

— [rarr) + TT {rcTzr) - tree = 0 , 
or oz 

the strain-displacement relations, 

ddr dr ddz 


d d 

^ (rCTzr) + ^ (rcjzz) = 0, 


1 f ddr 


^rr — 5 ^66 — 

or r 


^ZZ — 


Czr — 


ddz 

dr 


(3.1) 


(3.2) 


dz ’ 2\ dz 

and Hooke’s law, 

ECzr = (1 + , ECrr = Orr - V (Czz + OSe) , Eeee = CFeg - V {Orr + CTzz) , (3.3a) 

Eezz = CTzz -V {arr + CTgg) . (3.3&) 

The boundary conditions for the stress applied by the liquid at r = -I- dr are 

/ J 5 r. , {dUz dUr\ , 

arr(T = ri + dr) = -p + 2pL-^ , Gzr[r = Ti + dr) = -pL \ , (3.4a) 

and at r = To -I- dr the externally applied stress is 

arr{r = ro + dr) = —Pe , (Jzrir = ro + dr) = 0 ■ (3.4&) 
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Following elastic shell theory ( Mollmam]l 198 il l . we define stress resultants for the forces 


Uij and moments niij as 


riij = / (Tijdr, rriij = / a^rdr. 


(3.5) 


The curvature of the external shell dr^/dz will introduce only negligible terms, of order 
0(£i), in the stress resultants (13.51) . However, centre-body curvature drb/dz will have 
leading order effects in the fluidic analysis (§3.2). We multiply the axial momentum 
equation by r — and integrate over r to obtain a resultant form. We then differentiate 
with respect to z in order to relate to the radial momentum equation dSH). This yields 


dn. 




m. 


, da^r {ri) 

+ ri -5- = 0(e2). 


dz dz^ dz 

Substituting into the resultant form of the radial momentum equation, 


riP - ToPe + ri 


' d^m 


\ dz^ 


zz dcTzr {ri) 

- ri 


dz 


- nee = 0(^2). 


(3.6) 


(3.7) 


Order-of-magnitude analysis of the axial momentum equation in (13.11) yields a relation 
between characteristic values, 


nlz = drP*, (tIz = —£ 3, 
£2 

and order-of-magnitude analysis of (1X71) yields 

* 

— '£boP ) o'ee — • 

£2 


(3.8) 


(3.9) 


The shear and normal stress characteristic values, a*^. = p* and cr*^ = eip* are obtained 
from the fluidic problem (see §3.2). 

Taking note that m*„, = ri,Qn*^, normalizing (1X71) yields in leading order, 


Nee — P — Pe 0{eie, £ 2 ) 


(3.10) 


and normalizing Hooke’s law (13.31) accordin g to the ab ove characteristic values yields a 
reduced form of Love’s first approximation ( Love 1888fl . 


Ee2 Es2 

^99 ~ I ^zz ’ 

p* 


(3.11) 


for w hich orders of O { 62 , £ 3 ) are neglected and 'See = r^eelip*l£ 2 )- We apply ( Dugdale fc Rui j 


I 971 II the Kirchhoff hypothesis and describe the displacement field in terms of the radial 
dr and axial dz displacements of the midsection, denoted by overbars, 


- ddj- — 

dz=dz + {r-rm)-^. dr = dr , 


and thus we can represent the strain as a function of the deformation by 

dr 

„ eee = 

oz 


ddz , .d^dr 

£zz = — - {r -Tm)- 


(3.12) 


(3.13) 


dz'^ r 

Integrating (13.111) over r into resultant form and substituting Nee from (13.101) provides 
the leading order fluidic pressure to elastic deformations relation, 

fZ 


Dr = P-Pe, Dz-Dz{Z = 0) = -U {P-Pe)dZ, 


(3.14) 
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as well as a relation between the characteristic radial deformation and characteristic 


d*r 

= — = £3 < 1 • 


(3.15) 


pressure, 

P* 

E£2 no 

Equation (13.151) provides limitation on the maximal allowed pressure for which the as¬ 
sumption of small deformations is valid. 


3.2. The fluidic problem 

We assume an axisymmetric incompressible Newtonian fluid, governed by the momentum 
equations. 


dUr dUr dUr 


P 


duz 


duz 


du. 


dt ar dz 



1 d 

^ dUr\ 

d'^Ur 

Ur 


r dr 

y dr) 

dz^ 

J .2 

dp , 

ri d 

( duz^ 


+ pg 


r dr 


dz^\ 


and conservation of mass. 


Id, , duz 


The relevant boundary conditions are no-slip and no-penetration at r = dr, 

/ , , ^ ddr , , , XX ddz 

Ur{r = rb + ho+ dr) = , Uz(r = rb + ho + dr)) = , 

and at the centre-body boundary, r = Vb, 

Ur{r = Tfe) = 0 , Uz{r = Tb) = uo . 

Pressure at the inlet and outlet may be governed by 

p{z = 0) = po{t), p{z = 1) = pi{t), 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 


where po{t) and piff) are arbitrary functions of time. We normalize the coordinate s = 
r — rb{z), defined in the gap 0 ^ s ^ ho {z) + dr (z, t), by d*, 

S=-^, 0 ^ S XhHoiZ) + DriZ,T), (3.22) 

where Hq = ho/ho- Based on (I3.22p we can relate the normalized radial coordinate R to 

5, 

R = l + e3{S + XbSb{Z)), (3.23) 

where Sb (Z) = Sb (z) /si- Order-of-magnitude analysis of p.l8l) yields 


u 

u 


* 

r_ 

* 

z 



£1 , ul 


£id*p* 


(3.24) 


We require Uq ^ u* so that the viscous stresses resulting from centre-body motion scale as 
fiuo/d*. We focus on negligible gravity, G = pgl/p* ^ 1, and define a reduced Reynolds 
number. Re = pu*dr/ii, and Womersley number, o? = pd*'^/pt*. Applying transforma¬ 
tion (13.231) to (13.161) - (13.181) and employing (I3.24|) yields the following reduced system for 
the fluidic domain. 


dP 


d^Uz 

dS'^ 


O 


£3, £iRe, ef, ^b-^, G 


(3.25a) 
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O [a^,elRe,el ), 


(3.256) 


dUr , dSbdU, dU, 

dS ~ dZ dS ^ dZ ~ 


(3.25c) 


Equations (13.2511 are valid for Xh,Xb ~ 1 and Xh,Xb 1- t* represents the time-scale of 
the viscous-elastic interaction and is to be obtained by relation to the elastic problem. 


3.2.1. Nonlinear problem formulation 

The boundary conditions (13.191) . (I3.20p take the normalized form, 


Ur {S = XhHo + Dr) 


d*r dDr 
t*^~^ ’ 


U, {S = XhHo + Dr) 


dl dD, 

t*ul dT ’ 


and 


(3.26a) 


[/^(5' = 0) = 0, Uz{S = 0) = Uo, (3.266) 

with Uq = uoful and T = t/t*. We solve (I3.25al) . (13.2561) imposing conditions (I3.26al) . 
(|3.266l) . The resulting leading order axial speed reads 


S 

Ho + Dr)\ ’ 

(3.27) 

and is defined for XhHo + Dr > 0. Addressing mass conservation (I3.25cl) . we install (13.271) 
and integrate with respect to S across the film layer to produce a reduced Reynolds 
equation relating fluidic pressure to elastic deformations and curvature, 


U, 


1 dP 

2^ 


[S^-{XhHo + Dr) S] + 


S 


dl dD^ 


{XhHo + Dr) t*ul dT 


Uo 


1 - 


{Xh 


d*r dDr 

¥^r~^ 

1 4* 

^ 2 t*ul 


J_A 

12^ 


^{XhHoPDrf 


Uo 


Id dSb 

2^(A.i?o + i?.) + A,— 


{XhHo + Dr) 


d‘^D, 

dZdT 


^^{XhHo + Dr) + 2A{,—— 
dZ dZ 


dD, 


dT 


0 . 


(3.28) 


Equation (13.2811 is valid for XhHo + Dr > 0 and must be zeroed for values of {Z, T) that 
render XhHo + Dr ^ 0, since fluid pressure cannot contract the tube beneath the centre- 
body surface. We are now in a position to evaluate the characteristic coefficients and 
hence the time-scale t*. By dominant balance a consistent result is achieved by setting 
t* = dl/ul- The residual axial speed terms, coefficients of dD^/dT, are of order O {eo). 
Thus (13.2811 reduces to 


dDr 

~dT 


I d {dP 


12 dZ dZ 


;{XhHo+Dr)''\pUo 


19 dSb 

2^(A,.i7o + ^.) + Afc — 


0. (3.29) 


3.2.2. Linear problem formulation 

In the case of large Xh = 1 the characteristic gap should be defined by hg 

and not d*. This requires renormalization of the lubrication approximation attained in 
(|3.25all - (|3.25cll . The renormalized parameters associated with the linear problem, denoted 
by the tilde symbol, are given by In the case of large A/j = hg/d* 1 the characteristic 
gap should be defined by hg and not d* . This requires renormalization of the lubrication 
approximation attained in (13.25all - (13.25cl) . The renormalized parameters associated with 
the linear problem, denoted by the tilde symbol, are given by 

£i = Xh£i, £3 = XhEst Xb = Xf^^Xb, Re = XhRe, d = XhCx. 


(3.30) 
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We renormalize the local film coordinate and domain, 

S = Xl^S, ^Ho{Z) + \'^^Dr{Z,T), (3.31) 

and reapply boundary conditions (I3.26a|) - (I3.26&I) in S. The characteristic radial and axial 
velocities now scale as u*/u% ~ ei and u* = iih^p*/^. The reduced equations of motion 
for the linear case may be rewritten in tilde form, 

^ ^ + O ,h,eiRe,el,\b^,G^ , (3.32a) 


dP 






^-0(£3). 


(3.326) 


(3.32c) 


Following a similar procedure to that of §3.2.1 we derive a linearized reduced Reynolds 
equation analogous to (I3.29|) . 


dDr 

'W 


J__5 

12 ^ 1 dZ 


^^{Ho + X^^Dr)''\+Uo 


^ ^ {Ho + Xj;^Dr) +Xb'^^^ 


2dZ 


dZ 


0. (3.33) 


3.3. The fluidic-elastic problem 

We substitute (13.1411 into ()3.29p dropping the midsection overbars and derive an initial 
boundary value problem for the nonlinear viscous-elastic interaction (for the correspond¬ 
ing dimensional equation see equation lA II of Appendix [AT) . 


d{P-Pe) 1 d (dP 


:{XhHo + P-Pef\+UQ 




dT 12 dZ [dZ 
with corresponding centre-body surface condition rewritten in pressure form. 


0 , 
(3.34) 


P > Pe — XhHo . 


(3.35) 


Equations (I3.34|) - (l3.35p may be solved with an initial condition P{Z,0) = PflZ) and 
boundary conditions of type (13.211) in normalized form P(0, T) = Pq{T),P{1,T) = Pi(r). 
The characteristic time-scale t* of the nonlinear viscous-elastic interaction is evaluated 
via relations (13.151) and (13.241) . 


{E£2)^P 

,^*3^2 


(3.36) 


The dependence in characteristic pressure stems from the pressure dependent diffusion 
coefficient of (13.341) . 

As per the linear case, all elastic analysis results of §3.1 hold with substituted tilde 
parameters of (13.3011 where applicable. Therefore an analogous substitution of (13.141) into 
(j3.33ll will yield the corresponding linear initial boundary value problem for inlet pressure 
driven viscous-elastic perturbations of the external shell, 


d{P - Pe) 


dT 


12 dZ 


{^{Ho + X^\P-Pe)y^ 


+Uo 


1 d 


dSb 


2^ (^0 +A/(P-Pe)) + At,— 


dZ 


0 . 


(3.37) 


With corresponding initial and boundary conditions. The characteristic time-scale t* is 
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reevaluated for the linear case, 


t* = 




j,~2~ > (3-38) 

Ee{£^£2 

and shows exclusive dependence in solid-liquid material properties and the geometry of 
the configuration. A result which is c onsistent with the linear interaction time-scale of 
the full cross-section cylindrical shell (lElbaz fc Gatll2014ri . 


3.4. Non-linearity and neglected inertia terms 
Equating the non-linear time-scale, (13.361) and linear time-scale, (I3.38|l we derive a condi¬ 
tion on the characteristic pressure for which the system will exhibit nonlinear behavior, 

p* > EXh£2£3 = ^ . ( 3 . 39 ) 

Ti ri 


As the pre-wetting layer vanishes, hg ^ 0, any small perturbation in pressure will prop¬ 
agate nonlinearly, representing a state for which the problem is not linearizable. There is 
an upper limit on the characteristic pressure beyond which the assumption of negligible 
inertia (j3.25Q|l - (|3.256p will no longer be valid, representing the range of validity of the 
analysis. Substituting (|3.36p into (|3.25ajl . we demand that both and £iRe terms be 
negligible and derive the condition that 


p* < 



{E£2y^ 

P5 


(3.40) 


It will be shown in §5 that the transition from a nonlinear to a linear regime occurs quite 
rapidly in the intermediate Xh range. 


4. Closed-form solution of the peeling problem and derivation of 
propagation laws 

At the limit Xh ^ 0 while setting P^, = Uq = 0, (13.3411 describes cylindrical pressure- 
driven viscous peeling of the external shell from the centre-body and reduces to the 
porous medium equation (PME) in fourth power, 


dP d^P^ 
df ~ ~dZ^ 


(4.1) 


with adjusted time-scale, T = t/nt*, k = 48. We transform to self-similar variables. 


P = T-‘^f{v), V = ZT^ 


(4.2) 


and adopt Huppert’s renormalization (lHuDDertill982r) in light of the close mathematical 
analogy to a viscous gravity current propagating under a density gradient, 


E{0 = Vp^f{v) N = Vf^- 


(4.3) 


rjp describes the propagating peeling front for which the solution is supported in the 
region Z < Zp{T), prior to boundary interaction Zp{T) < 1. We install (14.21) . (14.31) into 
(14.11) which then reduces to the following ODE in F, 

(ic4)" + CE' + aF^O. (4.4) 


The eigenvalue a describes the inlet pressure boundary condition, 


p(o,r) = T-“, 


( 4 . 5 ) 
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and thus, 

i^(0) = r;;§, F{^-^1)=0. (4.6) 

We complete the formulation with integral mass conservation represented in non-dimensional 
form by 


P{Z,T)dZ = QT- 


Q = Vl 


F{Od(. 


(4.7) 


Both rjp and Q are functions of the eigenvalue a. Contrary to a viscous gravity current 
where the interface is induced by a given flux input, in which case the flux rate Q may be 
set constant for all a, in the current study the flux rate derives from the inlet pressure. 

Eq. (14.4p is solved analytically for a = 1/5 and is a particular case of the source solution 
of the PM E, known as the Z KB solution, as first obtained by IZeldovich fc Komnaneets 
( 1950( 1 and BarenblatH ( 1952h . The underlying boundary and initial conditions are 


P{Z,Q)=Q5{Z), (4.8) 

representing a sudden input of mass, Q, into the interface at T = 0, after which the inlet 
is sealed. The amplitude decay at the inlet, as the front spreads through the interface, is 
then given by The solution reads 


6 

where (s)^ = max(s,0) and T is Euler’s Gamma function. The velocity field correspond¬ 
ing to (14.911 can be attained via (13.271) and (13.25 ell . The dimensional solution correspond¬ 
ing to (14.91) is given by (IA2I) . For all other physical values of a (a < 1/5) (14.41) . along 
with conditions (1461) and (lO) . is solved numerically. In the range a ^ 0 the underlying 
initial condition is P{Z,Q) = 0. Fig. presents the self-similar pressure profile F{^) 
for a = 1/5 (|4.9I) and for a = 0, —1, —2, —3 corresponding to inlet signals of type (14.51) . 
Fig. [2Kb) depicts the flux rate Q, the front locus r]^, and the time it takes the front to 
reach the opposite boundary, denoted T;,, as a function of a. Fig. |2Kb) is plotted in the 
range a ^ 0 (the range 0 < a ^ 1/5 is irrelevant for comparison due to the change in 
boundary and initial conditions). Tf, represents the limit of validity of the self-similar 
analysis when adapted to a finite domain, it reaches a maxima around a = — 1 bel¬ 
low which it decreases moderately. Figs. n c) and (d) depict the resulting deformation 
regime, and Dr, respectively, for the case of constant boundary pressure (a = 0 ). 
The peeling front enters from the left {Z = 0) while the tube is set stationary at the 
right end Dr{l,T) = Dz{l,T) — 0. The deformations are attained via (13.141) . a strain 
ratio of £3 = 0.25 and a Poisson’s ratio of = 0.5 were used. As the front propagates the 
interface the tube shrinks from the free left end. The upper half of the radial deformation 
profile slightly steepens due to the deformation field. 

The propagation laws for an inlet signal of type (14.51) . according to self-similarity, are 
written in terms of the front location Zp, accumulated mass in the interface M, and time 
to boundary Tf,, 

Zp = r]p{a)T^~, M = Q{a)T^~, Tb = r]p{a)~^^ , (4-10) 



P{Z,T) =T- 


1 - Z^T-i 

40 


Q = 


-I + 


r(^ 


where T = t/nt* and k = 48. The corresponding dimensional propagation laws are 
given in (jA 31) . The spread rate of a disturbance in an unperturbed medium in the case 
of a full cross-secti on cylindrical shell containing a viscous liquid is of order 0(t^/^) 
( Elbaz fc Gat! 2014l l. and is independent of the inlet signal or the nature of the external 
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Figure 2: Self-similar solutions of cylindrical pressure driven viscous peeling, correspond¬ 
ing to Xh ^ 0,Pe = Uo = 0 and inlet pressure P{0,T) = T~°‘. Part (a): Pressure profile 
F(^) corresponding to eigenvalues: a = —1/5 (analytic) - dashed line, a = 0, —1, —2, —3 
- solid lines (obtained numerically). Part (b): Flux rate Q{a) - solid line, front locus 
rip{a) - dashed-dotted line, time to boundary Tt,{a) = tt/nt* - dashed line {Q, and Ti, 
obtained numerically). Parts (c,d): Constant boundary pressure viscous peeling (a = 0), 
tube fixed at Z = I, 0^(1,T) = D^(\,T) = 0. Part (c): Axial deformation D^{Z,T), 
progressive snapshots of AT = 1.75, between 0 ^ T ^ 15.75, shown in the natural 
time-scale T = t/t*. Part (d): Radial deformation Dr{Z,T) vs. undeformed axial coor¬ 
dinate - dashed black, vs. deformed axial coordinate - solid red, progressive snapshots of 
AT = 1.75, between 0 ^ T ^ 15.75, shown in the natural time-scale T = t/t*, ez = 0.25, 
V = 0.5. 
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forcing, this is a linear property which stems from the heat kernel and applies to the 
linear interaction model of the current study (j3.37|) . Thus, when comparing the spread 
rate of pressure driven viscous peeling with the propagation of a disturbance over a non- 
negligible fluid layer, we note that an impulse of type (14.81) is slower, than its 

linear counterpart for large enough time, and that otherwise faster rates can be achieved 
due to the dependence in inlet pressure. 


5. Effect of pre-wetting layer, self-similar analysis and numerical 
solution 


We consider a constant nonzero pre-wetting layer A^, > 0, i/o = 1, about a stationary 
cylindrical centre-body, with no external forcing Pg = 0. Eq. (I3.34p reduces to 


dT ^dZ 




0 , 


(5.1) 


with T = t/nt* and k = 48. We transform back to PME form, ^ = Xh + P, 

^ ^ ■ 

Focusing on the case of constant boundary pressure driven peeling. 


(5.2) 


P(0,T) = 1, P{Z,0) =0, a = 0, 


(5.3) 


the self-similar analysis of §4 may be extended to include the effect of the pre-wetting 

layer. We adopt the tilde symbol for all self-similar variables corresponding to a pre- 

2 ^ ^ 

wetting layer. Following §4, we substitute $ = fjpF{^), and rewrite the boundary and 
initial condition in shifted pressure $ in terms of the pre-wetted self-similar pressure 
profile F, 

F(0) = r7;i(l + A„), F(e>l)=7);§A^, (5.4) 

in similar fashion to (HU). Since we are interested in values of Xh for which P{Z,T) 
maintains its compact support we must add the condition that, 


dPjl) 

di 

Integral mass conservation Kn becomes, 



P{Z, T)dZ = 


qt-2 


Q = vl [ F{i)di- Xhfjp , 

Jo 


(5.5) 


(5.6) 


where Q and fj^ are functions of the base layer thickness Xh- The self-similar ODE 
(63), at a = 0 ) written in F{^) along with boundary conditions (15.41) . (15.51) and (15.6|) 
formulate an eigenvalue problem in Xh- It is solved numerically with initial guess from §4 
(a = 0). The solution is presented in Fig. [Sja) and starts from the outer pressure profile 
of the unpenetrated interface {Xh — 0) previously attained in §4. Subsequent profiles 
correspond to Xh = 0.2, 0.4, 0.6, 0.8 and 1. Fig. |3Kb) depicts the resulting flux rate Q, 
front locus 77 ^, and boundary time Tf,, in the nonlinear scale range 0 ^ A/j ^ 1. We 
thus derive propagation laws for a constant boundary pr^sure viscous-elastic peeling 
front when spreading over a base layer. The front location Zp, accumulated mass in the 
interface excluding the base layer M, and time to boundary Tj, are given by, 

Zp=fipiXh)TK M = Q{Xh)T^, fh = fif{Xh), 


(5.7) 









Figure 3: Self-similar solutions of cylindrical pressure driven viscous peeling over a pre¬ 
wetting layer for constant boundary pressure, corresponding to 0 ^ A/i ^ 1, Pe = 0 
and inlet pressure P{0,T) = 1. Part (a): Shifted pressure profile P(^) corresponding 
to eigenvalues: Xh = 0,0.2,0.4,0.6,0.8,1 (obtained numerically). Part (b): Flux rate 
Q{Xh) - solid black line, front locus f}p{\h) - dashed-dotted black line, time to boundary 
Tb{Xh) = tb/nt* - dashed red line in right alined scale (obtained numerically). 


where T = t/nt*. In the presence of a base layer, the unpenetrated interface front of 
§4 (a = 0) is accelerated by fjp{Xfi)/fjp(0), resulting in faster boundary time Tb and an 
added mass of Q{Xh)/Q{0). 

For nonzero Xh (EU is no longer parabolic degenerate for P = 0 and therefore does 
not exhibit the orthogonal fronts of the PME (14.11) . As Xh grows, a larger fj^ is required 
to converge to a constant value at the right boundary (^ = 1) causing the profiles to 
spread inward (see Fig.jSja)). For A/j > 1, P scales as. 


F{i;Xh) = 0 



A/i > 1 


(5.8) 


and ultimately vanishes as Xh oo. Thus, for A/i > 1 the self-similar profiles of F{^) 
still provide an exact solution for P{Z,T) albeit in a diminishing time frame T < Tb, 
representing the limit of the nonlinear interaction model when extended into the linear 
Xh range. We’de like to emphasize that while constant boundary pressure driven peeling 
propagates as 0(t^/^) irrespective of the base layer thickness, this is not the case for other 
types of inlet signals. For example, source type boundary conditions (to be examined 
below (15.91) 1 transition from propagation for A/j —0 to for Xh —>■ oo. 

We illustrate the effect of a pre-wetting layer on the propagation of a source type 
viscous-elastic front via numerical solution of EH) ( readjusted to the natural time-scale 
T = t/t*) along with appropriate boundary conditions. 


P{Z,0)=S{Z), 


dP{0,T) 


5P(i,r) 


(5.9) 


A no flux condition at the right boundary has been added to demonstrate the interaction 
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of the front with the opposing wall, deemed to start around T = KTb{a = 1/5), for 
small enough values of Xh- Propagation in time is shown in 3 progressive snapshots in 
Fig. |3Da),(b),(c) corresponding to T = 0.012,0.12,0.6, respectively, and for varying base 
layer thickness 0 ^ Xh ^ 3/2. The solution was obtained via a finite difference scheme 
and was validated on the basis of section §4. For A/j > 1 the solution is presented in 
the nonlinear time scale normalization (13.361) for comparison. Acceleration of the front 
as well as its spreading is observed as the base layer thickens. At the physical limit of 
large Xh the solution for source type boundary conditions (15.91) . representing pressure 
driven viscous-elastic perturbations of the external shell, tends to the behavior of the 
heat equation, A/i —>■ oo in (13.371) . which converges exponentially and for which the 
transfer of information is no longer confined to a compact support of the solution, which 
generally travels immediately to the whole interface as of T = 0. Thus the solution of a 
propagating front or disturbance over a pre-wetting layer can be viewed as the synthesis 
of two fundamental solutions, the ZKB profiles of (14.91) and the Gaussian profiles akin 
to (13.371) . As the base pressure falls below the crossover value of condition (13.391) the 
solution will increasingly exhibit linear behavior and vice versa. 


6. Dipole structures and their role in the transient response to 
moving external forces 


One of the fundamental solutions of the PME ()4.ip . as first attained bv IZeldovich fc Barenblatt 
(|l957t . is the dipole solution. It is constructed by taking the distributional spatial deriva¬ 
tive of the delta function as initial condition, P{Z,0) = d6{Z)ldZ, on the infinite line. 

The dipole solution h as also been studied in the context of viscous gravity currents 
( King fc WoodsI 2003h . analogous to the configuration presented in §4. However, the 
peeling formation of §4 cannot sustain negative pressure values due to violation of the 
surface condition (13.351) and therefore cannot accommodate the classical dipole solution. 

Let us examine governing equation p.34l) when forced by external pressure Pe■ A non- 
negligible fluid layer resides in the annular gap Xh = l,77o = 1 and the centre-body is 
stationary Uq = 0 and of zero curvature dSb/dZ = 0. Hence (I3.34|) reads. 


dP 1 d (dP 

~ nlz 


P-Pef \ ^ 

’ ' dT 


( 6 . 1 ) 


A structure similar to a dipole may occur when a high pressure region and a low pressure 
region form adjacent to one another while the diffusion coefficient (here {Xh + P — PeY) 
zeroes at the pivot (see Fig. [SJa)). Such conditions can be met at the transient phase of 
the response to an advancing spatially localized external force, 


P^Z.T) = 5{Z-1/2-ojT) . (6.2) 

Fig. [S](a) depicts the initial transient phase of the solution of equation (16.11) under these 
conditions, plotted between 0 ^ T ^ 0.012 in AT = 0.002 increments. The forcing (16.2|) 
is implemented as a unit Gaussian so as to scale locally as the film layer Xh^ its width 
0.1 and velocity w = tt, as depicted in Fig. [5Kb). Outside the dipole centre the base 
layer smooths out the pressure. As the external shock moves down the tube T > 0.012 
the dipole structure will collapse and the system will transition into a new pressure and 
deformation regime. 

While (16.21) is a specific case designed to create a dipole-like structure, dipoles are 
inherent to the initial stages of the response to externally moving forces, occurring in 
areas where the diffusion coefficient vanishes. We define a general external moving wave 
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Figure 4: Nonlinear viscous-elastic interaction over a pre-wetting layer - source solu¬ 
tion. Numerical solution corresponding to 0 ^ ^ 3/2, Hq = l,Pe = Uq = 0, 

0 ^ T ^ 0.6. Pressure plotted vs. axial coordinate. Boundary conditions of type (Ell). 
Xh = 0,1/3, 2/3,1, 3/2 from dark to light reds respectively. Part (a): solution at time 
T=0.012, Fig. Part (b): solution at time T=0.12, Part (c): solution at time T=0.6. 


force of the form. 


Pe{Z,T) = asm{kZ- ojT). (6.3) 

Fig.|6]demonstrates this transient dipole formation for the case of a = 1, fc = Ttt, co = IOtt, 
Pe plotted in Fig. [Oja) between 0 ^ 0.1. A no-flux boundary condition was employed 
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Z 



Z 


Figure 5: Initial transient phase dipole formation in response to an advancing spatially 
localized external force. Part (a): Fluid pressure dipole structure close-up, solution in 
AT = 0.002 increments between 0 ^ T ^ 0.012. Part (b): External pressure Pe = 6{Z — 
1/2 — uiT), implemented as a unit Gaussian, width 0.1, velocity to = tt, 0 ^ T ^ 0.012. 

at either end of the tube, 

3P(0.T) 3P(1.T) _ 

dZ dZ ' ' ' ' 

representative of a closed system. The solution (see Fig. (IKb)) is plotted between 0 ^ 
T ^ 0.1 in AT = 0.014 increments. The outermost dipoles are strongly influenced by the 
boundary walls. The asymmetry of the central dipoles results from the simultaneity of 
the transient buildup and diffusion processes. Thus, the dipole structure is a fundamental 
mechanism of the transient response to external forcing and is essential when considering 
a pulsating external environment. 


7. Transported ZKB profile, a building block for creating complex 
deformation patterns 

In this section we illustrate the use of the transport term of (I3.34|) dP/dZ ■ UqI2^ and 
the relation between the speed of pressure propagation and the gap (Xh), in order to 
obtain isolated moving deformations. The linear transport term of equation (13.341) may 
represent a relative motion of the centre-body, or alternatively, Uo{T) can be thought of 
as a slip velocity condition induced by an electric zeta potential or a similar mechanism. 
We set Pe = 0 and examine uniform geometry Hq = 1, dS^/dZ = 0, with a small 
pre-wetting layer A/j <C 1. Eq. (13.341) reduces to 


^ +P1 

dT 2 dZ l2dz\dZ^ ^ ’ 


0 . 


At the left boundary we employ an impulse sequence of the form, 


N 


p{q,t) = Y,a, 

n=0 


0 




-0 



(7.1) 


(7.2) 


where 0 is the Heaviside function. The signal is modulated in amplitude, and in 
width AT„ and can be sequenced non-uniformly by Tn- We set a zero initial condition 
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Figure 6: Initial transient phase dipole structure formation in response to an external 
moving pressure wave. Tube closed at both ends: dP{0,T)/dZ = dP{l,T)/dZ = 0. Part 

(a) : External pressure Pg = asm{kZ — ujT),a = l,k = = IOtt, 0 ^ T ^ 0.1. Part 

(b) : Fluid pressure in AT = 0.014 increments between 0 ^ T ^ 0.1. 


and a no flux condition at the right boundary, 

P(Z,0) = 0, ^^^^ = 0. (7.3) 

Under the above conditions a sequence of symmetric ZKB profiles will propagate through 
the interface. Assuming sufficiently small A/i, their compact support is maintained un¬ 
der transport, that is, each pulse is independent and generally does not interact with 
neighboring pulses. The rate of diffusion of the profiles decays strongly beneath a cer¬ 
tain amplitude (determined by (13.361) 1 and is significantly slower than the rate of their 
transport, that is, they maintain their form as they are transported. Using this approach, 
isolated fluid segments can be generated arbitrarily and transported along the interface 
independent of the pressure gradient and having limited or negligible communication 
with each other. In particular, the boundary sequence (ESI) can be modulated to form 
a moving wave signal, as illustrated in Fig. 0 via numerical solution of (I7.1I) - (I7.3I) . A 
pre-wetting layer of \h = 0.15 was used with a constant centre-body motion of Uq = 2. 
Fig. [7Kb) shows 3 progressive snapshots of the resulting radial deformation wave at times 
T = 0.9, 0.92, 0.94, as it propagates the interface. The corresponding fluid pressure con¬ 
tour is plotted in Fig. [TKa) and shows the uniformity of the sequence beyond an early 
generation time, T > 0.35. The obtained solution closely approximates the waveform, 

Dr{Z,T)^^[lp sin(iFZ - UT)], § = ^ 
where A = 0.1 and K = = IOtt. 
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Figure 7: Moving wave deformation pattern generated via a modulated impulse sequence 
A pre-wetting layer of \h = 0.15 is used upon constant centre-body motion of 
Uq = 2. Modulation laws in the interval 0 < T ^ 1; A„ = [0.37, 0.54, 0.54, 0.54, 0.54], = 

[0.003,0.02,0.0373,0.055,0.0732], Ar„ = [0.0005,0.0026,0.0036,0.0044,0.005], n = 1..5. 
Part (a): fluid pressure contour in space and time. Part (b): 3 progressive snapshots of 
the radial deformation wave at times T = 0.9,0.92,0.94, as it propagates the interface. 


8. Concluding remarks 

The analysis presented here assumes several small parameters, including geometric 
requirements such as £2 = {fo — 1 and dynamic physical requirements such as 

small deformations £3 = d*/ri,o <C and negligible inertia Rcr = ^ 1- 

While the geometric requirements are given and constant for a specific configuration, the 
physical requirements depend on, and limit the magnitude of, the characteristic driving 
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pressure. To illustrate the maximal allowable pressures and characteristic time-scales we 
examine several configurations with constant geometric ratios e = £2 = 0.1. For water 
(/r = 10“^[Po • s], p = 10^[Kg/m^]) as the liquid and rubber {E = 10®[Pa]) as the tube, 
we obtain p* = 10^[Pa], = 10‘*[s] in the case of rbo = l[m] and p* = 6.3 ■ 10"^[Po], = 

4[s] in the case of r^o = 10“^[to]. For silicon oil (/i = 10[Pa • s], p = 7.5 • 10'^[Kg/m^]) 
and rubber {E = 10®[Pa]), we obtain p* = 4.2 • 10^[Pa],P = 10®[s] in the case of 
rbo = \[m] and p* = 2.6 ■ 10®[Pa],P = 0.5[s] in the case of r^o = 10 “®[to]. Hence, a 
wide range of characteristic driving pressures and time scales can be achiev ed by varyini 
the p roperties of the configuration. Unlike linear viscous-elastic dynamics ( Elbaz fc Gai 
201^, the assumption of creeping flow in the current nonlinear problem can be achieved 
for any configuration, as long as p* and \h are sufficiently small. 

Future research may focus on pressure and deformation propagation in cases with 
varying pre-wetting layer thickness Hq(Z) and centre-body contour Sb{Z). While the 
governing equation for such configurations is presented in this work (13.3411 . the current 
study did not examine solutions for varying Sb{Z) and Ho{Z), these functions may allow 
for more elaborate control of the deformation field of the external shell. In addition, 
future research may include the effect of varying external shell thickness, which is similar 
to spatially varying surface tension, and may allow for transport of isolated deformations 
without external mechanisms such as presented in §7. 


Appendix A. Summary of results in dimensional form 

Nonlinear viscous elastic interaction - governing equation in fluid pressure correspond¬ 
ing to (IS31- 


d{p-Pe) rgp d \dpf Ee 2 ho 

dt 12p,{Es2f dz I dz V rbo ^ 


n d 

( E£2hQ \ 

Ee2 dsb 

2dz 

+ P Pe 

\ no y 

-1 

+ 


3 


0 . 


(Al) 


Closed-form solution of impulse driven viscous peeling according to self-similarity, corre¬ 
sponding to (|4.9I1 . 


piz,t) =p*^ 


«P(P£2)^ 


3p*5 ^ Kn{E£2f\ " 2 

40P ye® ) ^ 


(A 2) 


Inlet pressure driven viscous peeling - self-similar propagation laws corresponding to 

(ICTUl) . 




p*^eH 

KfJ.{E£2) 


Vf (a) ■ 


e \ M 


Q{a) 


, , , Kp.{E£2) , 

= p.3g2 "dM) 


(A 3a) 

(A3&) 
(A 3c) 


V is the interface volume. 
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